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Abstract 

We study the existence of Strange Nonchaotic Attractors (SNA) in the family of Harper maps, proving 
that they are typical but not robust in this family. Our approach is based on the theory of linear skew- 
products and the spectral theory of Schrodinger operators. 
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The study of the attractors of a dissipative dynamical system is a topic of great interest, be- 
cause these invariant sets trap the evolution of a large subset of the phase space and capture the 
asymptotic behavior. It has been known for a long time that attractors can be strange 1 1], i.e. geo- 
metrically complicated. The first examples of strange attractors were chaotic, i.e. with dependence 
sensitive on initial conditions [2]. In [3] were found strange attractors that are nonchaotic, and it 
has stimulated much numerical experimentation (see the review [4]) as well as rigorous analysis 



Our interest in this Letter is to show the existence and abundance of Strange Nonchaotic Attrac- 
tors (SNA for short) in the family of Harper maps. This is a family of ID quasi-periodically forced 
maps that many authors have suggested as a scenario in which SNA appear 
By SNA we mean here an invariant set that is a graph of a measurable and nowhere continuous 
function (it is Strange), that carries a quasi-periodic dynamics (it is Nonchaotic) and it attracts 
exponentially fast almost every orbit in phase space (it is an Attractor) 113 8fl . We prove that these 
SNA are typical but not robust in the family of Harper maps, in the sense that they exist for a 
positive measure Cantor set of the parameter space. 

In our analysis, we exploit the connections between (a) the dynamical properties of the Harper 
map (a ID quasi-periodically forced map); (b) the spectral properties of the Harper operator (an 
example of a quasi-periodic Schrodinger operator); (c) the geometrical properties of the Harper 
linear skew-product (a 2D quasi-periodically forced linear map). 

In recent years our knowledge of the spectral properties of the Harper operator, also known 
as the Almost Mathieu operator, and related quasi-periodic Schrodinger operators has advanced 
spectacularly. The progress made will be relevant to our approach. In particular, the connections 
between (b) and (c) have been successfully applied to the solution of the "Ten Martini Problem" 
full , on the Cantor structure of the spectrum of the Harper operator. 

The connection between (a) and (c) in similar models has been used to study the linearized 

o 

dynamics around invariant tori in quasi-periodic systems [15]. Specifically, the formation of SNA 
in this linearized dynamics is suggested to be a mechanism of breakdown of invariant tori llfll . 

A consequence of our approach is that neither arithmetic properties of the frequency of the 
quasi-periodic forcing nor localization properties of the spectrum of the Harper operator are crucial 
for the existence of SNA. 

The family of quasi-periodically forced dynamical systems under investigation in this Letter is 
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the family of Harper maps 



a — b cos (2n9 n ) — y n ' 
* ' 

f a ,b(.Vn,8n) 

n +i = 6 n + uj (mod 1), 



(1) 



where y 6 R = [— oo, +oo] and # G T = R/Z are the phase space variables, a, b are the 
parameters, and uj is the /regency (it is assumed to be irrational). 
Notice that a Harper map is a skew-product map 

Fa,b,u(y n , On) = (fa,b(y n , #n), d n + Uj), 

defining a dynamical system inRxT whose evolution from an initial condition (y , 9 ) is described 
by the M/i-power (y , 9 ) = (f { a N b \y , B Q ), 9 + Ncu), for N e Z. 

In a Harper map the parameter a is called the energy or the spectral parameter because after 
writing y n = x n -i/x n this family is equivalent to the family of Harper equations, which are 
second-order difference equations 

x n+ i + x n _i + b cos (2^(6*0 + nu)) x n = ax n . (2) 

These equations are physically relevant because they show up as eigenvalue equations of Harper 
operators (also known as Almost Mathieu operators), 

{H h ^ fio x) n = x n+ i + x n _i + b cos (2tc(9 + nu)) x n . (3) 

These are bounded and self-adjoint operators on Z 2 (Z) whose spectrum, that does not depend on 
6*o, describes the energy spectrum of an electron in a rectangular lattice subject to a perpendicular 
magnetic flux OH- 

The formulation of the second-order difference equation © as a first-order system is the Harper 
linear skew-product 



I Xn+1 ' 
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a — b cos (27i9 n ) — 1 
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{ X n -1 1 



(4) 



Vn+l M a , b (9 n ) 

n +i = 9 n + u (mod 1' 



3 



whose evolution is given by the Harper cocycle 



M a>b (9 N ^)...M a}b (9 ) ifiV>0, 
/ if iV = 0, (5) 

M-l(9 N )...M-l(9^) if JV<0. 

Note that O describes the evolution of the slope y n of vectors v n under the action of the linear 
skew-product ©. That is, © is the projectivization of @. 

To understand the dynamics of Harper linear skew-products it is important to know the growth 
properties of the solutions. The exponential growth is measured by the Lyapunov exponents which 
we now define. Given any nontrivial initial condition of the skew-product ©, (v , 9 ) with t> 7^ 0, 
the (forward) Lyapunov exponent for (vq, 9q) is the limit 



K,b,u (vo , 9 Q ) = lim — log I v N I , 



lim — Ioe 

N->+oc N 



(6) 



whenever the limit exists (in which case it is finite). If Vq = (x-i, xq) and yo = X-i/x then one 
can also define the (forward) Lyapunov exponent of the Harper map for the initial condition 

(Vo, 0o ) by 



lim — log 



lim — log 



dy N 



(yo,0o) 
dy 

™iil(2/oA 



where 



,W („, a \ _ d fci,b f a \ ^Ub f a 
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<L(vo,o ) 



dy 



dy 



An easy computation shows the relation 



(3a,b,iu(yO,9o) = —2\ a:btU} (Vo,9o). 



(7) 



Backward Lyapunov exponents are defined by replacing limAr_> +00 with limTv^-oo in the above 
formulation. 

Oseledec I19I showed that for almost every initial condition (i> , 9 ) the Lyapunov exponent 
exists and equals the averaged Lyapunov exponent 



lim — / lo 
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MflU9) 



which is never negative and exists by the Kingman subadditive ergodic theorem [20]. 

The case of the nonzero averaged Lyapunov exponent, X a ,b,oj > 0, which we call hyperbolic, is 
important for our purposes. In this case, there exists a full measure set 9 C T such that for every 
8 G one has a splitting 



W s {8) © W u (8) (8) 



characterized by 



veW s (8)\{0} lim — log \M^ N \d)v\ = -A (9) 

and 

v G W u {8) \ {0} <S> lim — log \M {N) (6)v\ = +X (10) 

iV— >±oo A^ 

W s (#) and W / "(6 ) ) are the stable and unstable subspaces at 8, respectively. The elements of the set 
are referred to as the Lyapunov regular points. 

In the phase space M 2 x T, one can form the product sets W s and W u whose elements are 
pairs (v, 8) with v G W s (8) or W u (8) respectively (whenever these subspaces are defined). These 
are the stable and unstable subbundles. According to Oseledec lfl9Tl the ^-dependence of the 
decomposition is measurable but not necessarily continuous. 

When the splitting ® is defined for all 8, that is = T (hence ^-dependence of the subbundles 
is continuous B39J) the linear skew -product is said to be uniformly hyperbolic. Otherwise it is said 
to be nonuniformly hyperbolic. 

In the Harper map, we can determine whether or not hyperbolicity is uniform by looking at the 
spectral problem of ©. Indeed, an energy a is in the spectrum of the Harper operator © if, and 
only if, the corresponding linear skew-product © is not uniformly hyperbolic J^l^ . We will 
use an implication of this result: if a is in the spectrum of the Harper operator and the averaged 
Lyapunov exponent is nonzero at a, then the linear skew-product is nonuniformly hyperbolic. Let 
us now see that in this case the corresponding Harper map has a SNA. 

The above concepts of hyperbolicity can be translated to the dynamics of the Harper map (HJ) 
(which reflects how the linear skew-product © changes directions of vectors). Recall that, in the 
hyperbolic case, (3 < 0, there exist two invariant subbundles W s and W u , for 8 E 0, which are 
measurable as a function of 8 and satisfy © and (TTOb . We define y s (8) and y u (8) as the slopes 
of the subbundles: they are the only elements of R such that (1, y s (8)) T and (1, y u (8)) T belong to 
W s (8) and W u {8) respectively. The product sets 

Y s = {(y s (8), 8), 8eQ} and Y u = {(y u (8), 8), 8 G 0} 



are invariant under the Harper map and have quasi-periodic dynamics; thus Y s and Y u are non- 
chaotic invariant sets. 

Still in the hyperbolic case, the decomposition of M 2 into direct sum of W s {9) and W u (9), 
9 G O, implies that every pair (v, 9) of this skew-product @ other than the stable subbundle, is 
attracted to the unstable subbundle and grows exponentially in norm. Looking at directions (which 
is what the Harper map retains), forward orbits with initial condition (y, 9) (other than (y s (9), 9)) 
are exponentially attracted to Y u , that is 

1 



lim — lo£ 

iV->+oo N 



p = -2\<0, (11) 



while backward orbits (other than (y u (9), 9)) are exponentially attracted to Y s . Thus Y u is a non- 
chaotic attractor for the Harper map: for almost every initial condition, orbits are exponentially 
attracted to it. Similarly Y s is a nonchaotic repellor. 

Let us now relate the uniformity of hyperbolicity in the skew-product to the regularity of these 
nonchaotic attractors. If the skew-product © is uniformly hyperbolic then the invariant subbun- 
dles W s , W u are defined in all T and are continuous, and so are their projectivizations y u and 
y s . 

In contrast, if the skew-product is nonuniformly hyperbolic, then the invariant subbundles are 
measurable but not continuous and their projectivizations y u and y s are measurable but not contin- 
uous functions of 9. Moreover, discontinuities are propagated by the quasi-periodic dynamics and 
the invariance property of the attractor: if y u is discontinuous at a single 9 then the same happens 
for 9 N = 6>o + Nuo for all N G Z, so that the function y u is nowhere continuous. The same result 
happens for y s . 

In summary, in the nonuniformly hyperbolic case, we will say that Y u is a Strange Nonchaotic 
Attractor (SNA) of the Harper map because the following properties are satisfied: 

(i) Y u is the graph of a measurable function of 9, y u , which is nowhere continuous (Y u is 
Strange); 

(ii) Y u is an invariant set of the Harper map with quasi-periodic dynamics (Y u is Nonchaotic); 

(iii) Almost every orbit in phase space is attracted to Y u at exponential rate (Y u is an Attractor). 

Qi 

The existence of nonuniformly hyperbolic linear skew-products was already shown by Herman 
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1 2311 . who proved that 

A a)6ia ,>max(0,log^l ] (12) 



as long as uj is irrational. Moreover, Bourgain & Jitomirskaya 12411 prove that the equality in (fT2l 
holds if, and only if, a is in the spectrum of the Almost Mathieu operator. 

Thus, for \b\ > 2 and uj irrational, a Harper map has a SNA if, and only if, a belongs to the 
spectrum. Since the measure of the spectrum is given by the formula |4 — 2|6|| W25I1 these SNA 
are persistent in measure in the family of Harper maps. They are not, however, persistent in open 



sets, since the spectrum is a Cantor set 114 
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2711 . and therefore any SNA in a Harper map can 



become a regular attractor by means of an arbitrarily small perturbation. As an example, using the 
symmetry of the spectrum, the Harper map has a SNA for any |6| > 2, uj irrational and a = (this 
value always belongs to the spectrum). 

As an illustration of the above rigorous results, we performed several numerical computations. 
In the following, we chose the irrational frequency uj = ^ J41I1 ,6 = 3, and we considered a as a 
moving parameter. The averaged Lyapunov exponent as a function of a is displayed in Figure [T] 
Notice that STUi implies that A a > ^jp, so that the Harper cocycle is hyperbolic for all the values 
a. Moreover, the equality holds only if the cocycle is nonuniformly hyperbolic. As a result, the 
values of a for which Y u is a SNA of the Harper map correspond to the "flat pieces" of the graph 
in Figure [U which lie in a Cantor set of measure 2. The "bumps" appear in gaps of the spectrum, 
that is energies a in the resolvent set, for which Y u is an invariant attracting continuous curve. 
Hence, gaps are labelled by a topological index which is the number of turns of Y u on R 112 811 . 

We also selected several values of a and computed the attractor Y u and the repellor Y s of the 
Harper map for several values of a. The results are displayed in Figures |2] and |3] Since y E R 
represents a slope, in the pictures we display the angle cp = tan -1 y E [—ir/2, ir/2) corresponding 
to such a slope. In this representation we identify ip = ir/2 with ip = —n/2, because a line of slope 
+00 is a line with slope —00. The objects shown are graphs of functions y u (9), y s (9), respectively, 
of the slope y E R with respect to 9. 

In this Letter we have shown the existence of SNA in Harper maps with uj irrational, \b\ > 2 and 
a in the spectrum without resorting to the possible localization properties of the spectral problem. 
Recall that a is a point eigenvalue of a Harper operator H b ^ )9o if the corresponding eigenvalue 
(Harper) equation has a nontrivial localized solution ip = (i> n )nez which is square integrable or 
even decays exponentially with |n|. If the set of localized eigenvectors of a Harper operator forms 
a complete orthogonal basis of Z 2 (Z) then the spectrum is pure-point. 

1 




Figure 1 : The averaged Lyapunov exponent as a function of a. 

In previous work on the existence of SNA in Harper maps, localization was seen a s a j ustifica- 
tion for the strangeness of SNA, in the regime of nonzero Lyapunov exponents JsH^llJ 
As we have seen, we do not use localization to prove the existence of SNA. Besides, localization 
may not hold in all the Harper maps studied here. Indeed, an energy a in the spectrum of a Harper 
operator with nonzero Lyapunov exponent (for which the Harper map has an SNA) may not be an 
eigenvalue of the operator. Indeed, if cu is not Diophantine, localization may only hold for |6| > 2 
large enough (depending on uj) II27I1 . Even in the Diophantine case, the spectrum also contains a 
residual set of energies which are not point eigenvalues ||3i , 32 ]. Thus, there are SNA in the family 



of Harper maps without localization for the corresponding operator or arithmetic properties on ui. 
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a= 0.3, b= 3 




a= 0.4, b= 3 



9 o 




-nil 



a= 0.5, b= 3 





Figure 2: The attractor Y u and the repellor Y s of the Harper map for a = 0.3, 0.5 (SNA) and a = 0.4, 0.6 
(continuous invariant curves). Notice that a = 0.4 and a = 0.6 correspond to different gaps which are 
labelled by the indices of the continuous curves 8 and 5, respectively. 
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